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Abstract 



We propose a new boundary-driven phase transition associated with vor- 
tex nucleation in mesoscopic superconductors (of size of the order of, or larger 
than, the penetration depth). We derive the rescaling equations and we show 
that boundary effects associated with vortex nucleation lowers the conven- 
tional transition temperature in mesoscopic superconductors by an amount 
which is a function of the size of the superconductor. This result explains 
recent experiments in small superconductors where it was found that the 
transition temperature depends on the size of the system and is lower than 
the critical Berezinski-Kosterlitz-Thouless temperature. 
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Phase transitions in two dimensions (2D) has been the subject of long-standing interest. 
The main reason is that the Berezinski-Kosterlitz-Thouless (BKT) transition fl|-|3| must 
arise in 2D systems like, for example, superfluids and quasi-2D superconductors. Berezinski 

was the first to show that topological defects such as vortices may play a significant 
role in phase transitions. The fact that the energy required to create a vortex depends 
on the size of the system was the main obstacle to this fundamental discovery. Berezinski, 
however, recognized that the creation of a vortex-antivortex (V-A) pair is energetically more 
favorable, because the pair energy depends on the separation distance of the pair only. The 
mechanism of the BKT transition is: the creation of V-A pairs decreases the superfluid 
density, which in turn decreases the binding energy of the V-A pair. As a result, it becomes 
easier for more pairs to be created and again this reduces the superfluid density further. 
This renormalization process continues until the screening effect is large enough for vortices 
and antivortices to nucleate freely and spontaneously. Since single vortices destroy the phase 
correlation needed for 2D superfluidity, the superfluidity is destroyed and there is a jump 
discontinuity in the superfluid density at the transition. 

In three dimensions (3D), the creation of a V-A pair requires an enormous amount of 
energy depending on the size of the system and it was believed that a generalization of 
the BKT transition to 3D systems was not possible. However, Williams M and Shenoy |5j 
independently recognized that pairs of vortex loops may play the same role as V-A pairs in 
3D. With these ideas, Williams [|J constructed a model of the superfluid transition in 4 He. 
An analogous theory was developed by Shenoy || to describe the phase transition of the 3D 
XY model. More recently, Kusmartsev || proposed a mechanism for vortex nucleation in a 
flow of rotating superfluid 4 He, based on a mechanism similar to the BKT transition, where, 
however, half-vortex rings (HVR's) play a role similar to V-A pairs. The HVR's penetrate 
the barrier with the help of critical fluctuations via the creation of an HVR "plasma" . This 
leads to a topological phase transition where the barrier for vortex nucleation disappears 
and spontaneous threshold vortex generation starts. Kusmartsev derived the scaling laws in 
the critical region and introduced and estimated the scaling relation and scaling exponent 
for the critical velocity. 

The properties of a superconductor are expected to change radically when its size be- 
comes comparable to that of the Cooper pairs since the creation energy of a vortex is then 
of the same order as that of a V-A pair. Recent progress in microfabrication techniques has 
made it possible to study experimentally (mesoscopic and nanoscopic) superconducting sam- 
ples of micrometer and nanometer dimensions [|7|-[5|. This has led to the discovery of new 
phenomena in the superconductivity of mesoscopic systems, particularly the discovery of 
the paramagnetic Meissner effect [^-|T2[ • Moreover, while the value of the Ginzburg-Landau 
parameter A/£ (ratio of the coherence length £ to the penetration length A) is sufficient to 
determine the type of bulk superconductors, both experiments |||| and theoretical work 
T3HT7| on the magnetization of mesoscopic discs have shown that both the type and order 



of the transitions between the superconducting and the normal states depend on the size of 
the disc. 

In this Letter we propose a boundary-driven phase transition associated with vortex 
nucleation in mesoscopic superconductors not previously reported in the literature. The 
nucleation of free vortices from the boundary drives the system to a topological phase tran- 
sition with a lower critical temperature than the conventional or BKT critical temperature 



2 



of the system. Nucleation of single vortices is prevented by their attraction to the bound- 
ary — vortices induced in the system, just like 2D electric charges, are attracted to their 
mirror antivortices (vortices of opposite polarity). A single "free" vortex can penetrate into 
the system only by overcoming the vortex-image antivortex (V-IA) interaction or, in other 
words, a surface or Bean-Livingston barrier. However, creation of other V-IA pairs close to 
the boundary may renormalize the "Coulomb" attraction to the boundary: the V-IA plasma 
screens the attraction of the vortex to the boundary just like in the BKT transition where 
the creation of V-A pairs screens an effective interaction between a vortex and an antivortex. 
Eventually this leads to the creation of "free" vortices. These penetrate into the system and 
the order parameter associated with superconductivity is destroyed. 

To illustrate this effect, let us consider a circular flat superconducting disc of radius 
R = 0(A) and thickness d <C R and a vortex of vorticity k = h/m* (m* is the mass of the 
Cooper pairs) at a distance r < R from the center of the disc. Then it is straightforward to 
show that the image vortex (of vorticity —k) is a distance r' = R 2 /r from the center of the 
disc on a straight line joining the center to the vortex at r (see, for example, Ref. fl8j). The 
interaction energy Uq of the V-IA pair depends logarithmically on the separation r' — r, 

R 2 - r 2 
U = 2q 2 ln + E C , 

rr c 

where E c is the potential energy associated with the core of the vortex and r c is the effective 
core radius. By analogy with the 2D Coulomb gas, 

Ps\ 1/2 \l/2 ( h 



is the effective vortex charge. p s = p 2 ® = pf^d is the 2D superfluid density. U${r) is the 
energy with which the vortex charge at r is attracted to the boundary (surface) of the 
superconducting material. 

At low temperatures, near T = K, it is unlikely that more than only a few vortices 
will be present. However at higher temperatures, there are likely to be many more vortex 
excitations, including some located in the space between r and R. These have an attenuating 
effect on, and screen, the interaction Uq(t). Following Kosterlitz and Thouless, [pffl and 
Williams J| and Shenoy ||, we take into account this screening effect by introducing a 
scale-dependent dielectric constant 

e(r) = l + 47r X (r). (1) 

The effective susceptibility x( r ) = J R ot(r')dn(r'), where a(r) = q 2 (R — r) 2 /2kBT is the 
polarizability and n(r) is the number density of vortices. It is straightforward to show that 
dn(r) = 2irrdrexp(—U(r)/kBT)/r^. U(r) is the screened interaction, 

, r<"'-'V* dr' 

Introducing the dimensionless superfluid density K = q 2 /(irksT) and the renormalized 
density K r = K/e(r), Eq. (1) takes the form 
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k: 



K~ 



/ dr(R — r) 2 rexp —2-kK I 

Jr c Jr, 



~{R 2 -r 2 )/r fal 

e(r') r' 



where yo = exp(— E c /kBT). This derivation implicitly assumes a rather low density of 
vortices. This is evident, for example, in the fact that we have used the unrenormalized 
charge q instead of q r = q/e(r) to determine the polarizability. Next we also neglect the 
correction term in the V-IA interaction energy. Although the principal result is not changed, 
these approximations are necessary to prevent the equations from being intractable and lead 
to 



k: 



47r 3 ?/o 



R 



dr(R — r) 2 r exp 



-2nKln 



R 2 



rr r 



(2) 



To make Eq. (2) self-consistent, one needs to replace K in the exponential by the renormal- 
ized density K r . However, at low temperatures, the integral is small, and Eq. (2) is the first 
two terms in the expansion of K~ l . 

At temperatures near the phase transition where the superfluid density tends to zero, the 
perturbation series is not valid. In this regime, we use the vortex-core rescaling technique 
proposed by Jose et al. [^UJ (see also the review by Wallace |2U[]). The procedure is to split 
the range of integration [r c , R] into two parts, [r c , br c ] and [br c ,R], with b — 1 == Info <C 1, 
and only evaluate the non-singular contribution of small r. Rescaling r — > rb in the second 
integral to restore the original cut-off r c , we find a perturbative expansion for K~ l : 



-K: 



+ 



l 

fo 

4vr 3 y 



47T 



K- L + — r 2 c (R - r c y y exp -2nK In 



R 2 



R/b 



dr 



R 



— r 



rb 3 exp (—2nK In b) exp 



Info 



-2vrinn 



R 2 /b 2 



rr r 



(3) 



We now require that Eq. (3) has the same functional form as Eq. (2). This is achieved by 
introducing new variables at the increased scale, 



K 



i-i 



K~ l + 4tt : 



-yln fo 



y' = fo 3 y exp (-2nK Info), 



(4a) 
(4b) 



together with R' = b~ 1 R and K r = bK r , where we have introduced the fugacity y = 
yo exp{— 2itK ln[(i? 2 — r 2 )/rr c ]}. In deriving the above transformation, we have only retained 
terms of 0(y). It is convenient to build up a large increase in the core radius r c by successive 
repetition of this transformation. In this way, one arrives at differential renormalization 
group equations for the effective couplings K\ and yf. 



dKi 
dl 



K x - 4vr 



, (^-r c ) 2 2 

~2 K l Vh 



d i = $-2.K l)yh 



(5a) 
(5b) 



with the definition dl = Info. The scaled radius Ri of the disc satisfies dRi/dl = —Ri 
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The fixed point of the rescaling equations is defined by 
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(3 - 2nKi) yi = and K t - An 3 ^ ^ K? Vl = 0, 
which have the nontrivial solution 

The critical point (K*,y*) separates the two phases of the system: the first is the low 
temperature phase, characterized by growing superfluid (renormalized) density K\ m K e l 
(with K being the initial value of K t at scale size r c ) and vanishing fugacity y t e~ l ^°, 
£o = (tt-^o) -1 - The second one is the high temperature phase, characterized by exponentially 
growing fugacity: for infinite temperatures, K\ — gives y% = yoe 31 , i.e., vortices proliferate. 

The phase transition between these two regimes is easy to understand: At low temper- 
atures, T = 0+K, there are only a few vortices in the system (small fugacity). These are 
attracted to the boundary and cannot nucleate. The coherence length £ ~ O(r ). As the 
temperature increases, there is a growing number of vortex excitations and these screen the 
attraction to the boundary. The superfluid density decreases, the scaled radius Ri decreases 
and the coherence length increases (as ~ e l ). As the temperature increases further, there 
comes a point at which the screening is large enough for vortices to nucleate freely. R[ ~ £, 
the scaling stops and a phase transition occurs. 

To find the behaviour of the scaling near the critical point, we rewrite the rescaling 
equations (5) in terms of scaled deviations from the fixed point: we expand K\ and yi 
around K* and y* as Ki = K*(l + k') and yi = y*(l + y'). The scaling equations (5) then 
become, to first order in k' and y', 

i(;)-(=i-i)C) 

where we have dropped the primes. Expanding k, y and R in eigenstates A±(l) = A±e x±l 
of the fixed-point stability matrix above, the eigenvalues are A + = (vT3 — l)/2 and A_ = 
— (VT3 + l)/2. These define the relevant and irrelevant axes in the K { — y t plane. We 
assume, following existing procedure, that the relevant scaling field A + is the temperature 
axis, A + Ki A\e\, where e = (1 — T/T c ) is the deviation of the temperature T from the 
transition (superconducting or Berezinski-Kosterlitz-Thouless temperature) T c , and A is a 
constant. 

Then the rescaling law for the free energy F per unit area implies 

Z(K , y ,R ) = e-^-^Z^uyuRi) 

= e-W-Wz (A\e\e x + l ,A_e x - l ,R e- 1 ) , 

where Z is the partition function and R is the (unsealed) radius of the superconducting 
disc. The scaling stops when R t = R e~ l reaches a critical radius R c = £ = r c e l . This 
happens when 

l = l c = \\n^. (9) 



(8) 
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Setting I = l c in Eq. (8), the partition function is well defined only if 



f(= r c e c ) = r c |e 



1/A+ 



(10) 



Combining Eqs. (9) and (10) gives |e| 2 = (r c /R Q ) x+ , whence 




) 



A+/2 



(11) 



This result implies that the superconductivity breaks down at a temperature which is lower 
than the conventional superconducting or BKT critical temperature T c . The mechanism of 
vortex nucleation into the disc we have described leads to a depression of T c by an amount 
which varies inversely with the radius of the superconducting disc 



This depression in the critical superconducting/BKT temperature has been observed in the 
recent experiments of Geim et al. |2l| on the magnetization of mesoscopic superconducting 
discs of various radii, typically with d ~ 0.1 /^m and Rq ~ 1-10 /im. For these (aluminium) 
discs, r c (equal to the coherence length £ & t zero temperature) can be estimated as r c ~ 
£o ~ 0.18 /im. 

The phenomenon described above is very general. Obviously it may arise in any small 
system, like superfluid droplets, quantum dots in a superconducting state and so on. It 
may arise not only in superfluid or superconducting systems, but also in other condensed 
states, such as the magnetic state. In any case, topological defects originating from the 
surface/boundary (like vortices in superffuids and superconductors) have the potential to 
destroy a condensed state provided that the system has a small size and, in general, the 
critical temperature must decrease with the size of the system. 

We have described a mechanism of vortex nucleation in superconducting systems of size 
comparable to the characteristic size (coherence length) of the quasiparticles of the systems, 
the Cooper pairs. A vortex created in a small superconductor is strongly attracted to the 
boundary of the system and hence cannot nucleate. The attraction is due to the Coulomb 
attraction of the "vortex charge" to the "vortex charge" of the image antivortex. However, 
a fluctuating creation of several of these topological defects in the bulk of the supercon- 
ductor, in the space between the vortex and the boundary, screens the V-IA attraction by 
renormalizing the superfluid density. This further improves the condition for the creation of 
more of such fluctuations. The screening effect of the plasma of vortex fluctuations continues 
until vortices nucleate freely. This leads to a phase transition at which the order parameter 
associated with the superconductivity of the system is destroyed. 

We have used a real-space renormalization group method to derive the scaling laws in the 
vicinity of this phase transition, and shown that the phase transition occurs at a temperature 
which is lower than the superconducting or BKT temperature T c . The amount by which T c 
is lowered is equal to T c (r c /_R ) a (with a > 0), showing that the depression in the transition 
temperature increases with decreasing disc radius R$, in agreement with recent experiments 
pl| on mesoscopic superconductors. 




(12) 
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The properties of a mesoscopic superconductor are remarkably different from those of 
macroscopic or even microscopic superconductors, being dependent on the size of the system. 
In this Letter, we have predicted a new phase transition that occurs in such confined systems 
and driven by the boundary of the systems. More precisely, this phase transition, although 
similar to the BKT transition, is driven by the creation of a plasma between the vortex 
and the boundary which screens the Coulomb attraction of the vortex to the boundary. In 
this respect, the driving mechanism for the phase transition is analogous to the mechanism 
proposed by Kusmartsev || for the nucleation of vortices in a flow of rotating superfluid 
4 He. 

We would like to thank A.K. Geim for discussing his recent experimental results with us 
prior to publication. We would also like to acknowledge financial support from the EPSRC. 
FVK is an associate fellow of the L.D. Landau Institute for Theoretical Physics, Moscow, 
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